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Abstract 

To investigate the properties of c = 1 matter coupled to 2d-gravity 
we have performed large-scale simulations of two copies of the Ising 
Model on a dynamical lattice. We measure spin susceptibility and per- 
colation critical exponents using finite-size scaling. We show explic- 
itly how logarithmic corrections are needed for a proper comparison 
with theoretical exponents. We also exhibit correlations, mediated 
by gravity, between the energy and magnetic properties of the two 
Ising species. The prospects for extending this work beyond c = 1 are 
addressed. 
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1 Introduction 



There is at present considerable analytic understanding of how conformal 
matter with central charge (c) less than or equal to one couples to two- 
dimensional gravity Jl]. These c < 1 models are relevant both as non-critical 
string theories and as novel statistical mechanical systems describing matter 
on a dynamical substrate lattice. The relation of these models to quantum 
Liouville theory allows one to compute the shift in the scaling dimen- 
sion of physical operators due to the coupling to 2d-gravity. This shift is 
determined solely by the scaling dimension of the operator in the absence of 
gravity (i.e on a fixed lattice) and the central charge of the theory. Given the 
dressed scaling dimensions one can compute critical exponents of thermody- 
namic observables related to correlation functions. Through the mapping of 
these systems onto matrix models it is even possible to incorporate topology 
change by performing the entire sum over topologies in the so-called double- 
scaling limit. This beautiful state of affairs falls apart for central charge 
exceeding oneQ. In this case the methods of Liouville theory appear to fail - 
in particular they predict complex critical exponents. The vanishing of the 
mass of the dressed identity operator (the tachyon) at c = 1 suggests the 
onset of an instability in the worldsheet geometry. Based on this, it is widely 
believed that c > 1 no longer describes a continuum theory of surfaces 
Since no models have been solved in this regime, it would be valuable to 
apply numerical techniques to see if c > 1 models appear to be qualitatively 
different than those with c < 1. 

A simple class of models to test our understanding of these issues is 
provided by multiple copies of the Ising model on a dynamical lattice. Each 
Ising model has c = |. A model with n individual copies then has c = n/2. 
For a single Ising model (n = 1) on a dynamical lattice the critical exponents 
may be computed analytically via a mapping of the model onto a particular 
2-matrix model (with one matrix representing spin up and the other spin 
down) |5], ||. The model has a third order rather than a second order phase 
transition. Without the coupling to gravity the partition function of the 
n-Ising model would simply be the n-th power of the single Ising model 
partition function, since the spins of different copies are independent. The 

1 More precisely when the quantity c — 24 A > 1, where A is the conformal weight of 
the lowest- weight state in the theory ||. 
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interaction with gravity, however, induces an effective interaction between 
spins of different species (copies). Spins of different species effectively interact 
through the dynamical lattice. Brezin and Hikami J7|, |8], || have studied an 
equivalent 2 n -matrix model realization of this system in perturbation theory 
in the cosmological constant, without any apparent anomaly setting in at 
c = 1. Similarly Monte-Carlo studies of multiple-Ising and multiple-Potts 
models on dynamically-triangulated random surfaces (DTRS) do not uncover 
dramatic changes as c passes through one [[H], [11], [12|] . 

On the numerical side it is important that we understand the transition 
case c = 1 before plunging into the regime c > 1. This is the motivation 
for the work presented in this paper on large-scale DTRS simulations of the 
n = 1 and n = 2 Ising model coupled to 2d-gravity. In the n = 2 case the 
formalism of KPZ allows a computation of the relevant critical exponents, 
which we can consequently compare to the numerical results extracted from 
finite-size scaling and direct fits. The exact solution of the single-Ising model 
also provides a direct check that the discrete DTRS algorithm is reproduc- 
ing continuum behavior for the lattice sizes simulated and that the numerical 
analysis methods employed are adequate. The n — 2 (c = 1) case is consider- 
ably more complicated than the n = 1 case because of logarithmic violations 
of scaling. Here the extensive work on c = 1 matrix models provides an 
essential clue ||13|| . It allows us to identify the appropriate scaling variable 
that replaces the cosmological constant. With this in hand we are able to 
compare the results of our large-scale DTRS Monte-Carlo simulations with 
the predictions of KPZ. 

Similar issues are also addressed in a companion paper |14| in which per- 
colation coupled to gravity with c = 0, 1/2, 1 and c > 1 matter is examined 
in much detail. The nature of finite size effects in simulations of two dimen- 
sional gravity plays a central role in that paper, as it does here. We shall 
show how the percolation results for Gaussian matter in [14] are consistent 
with those obtained in the two species Ising model discussed here. 

The outline of the rest of the paper is as follows. In section 2 we give a 
theoretical discussion of c = 1 conformal matter coupled to 2d-gravity, with a 
derivation of the critical exponents we will be comparing with our numerical 
simulations and a discussion of logarithmic violations of scaling. In section 3 
we present our numerical methods and results including a comparison with 
theoretical expectations. Finally, in section 4, we conclude with a discus- 
sion of the origin of logarithmic corrections to scaling and the prospects for 
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extending this work to the regime c > 1 that this implies. 

2 Theoretical Predictions 

We shall consider a model in which Ising spins are attached to the vertices 
of triangulations. The triangulations are characterized by their adjacency 
matrix which equals 1 if % and j are neighbors and vanishes otherwise. 
Cij is the discrete analogue of the worldsheet metric g^. We shall restrict 
ourselves to the set of triangulations with N vertices Tn containing only 
loops of length 3 or greater and vertices of coordination number of at least 
3. We simulate a theory determined by the partition function 

n s N 

E ex P (-/3]r £ Cv{T)o?of); (1) 

TGTjv CTi=±l a=l i,j=l 

a labels the spin species. In this paper, we address the cases n s — 1 and 
2. Most of the relevant theoretical calculations are performed in the grand- 
canonical ensemble, with the partition function 

Z(ji) = Z N exp(-fiN) (2) 

N=l 

dependent on /i, the cosmological constant. 

Our primary observable will be the spin-susceptibility, which we express 
as the integrated spin-spin correlation function 

^ = ^(EE« (3) 

ll s iv a= i ij 

The integrated spin-spin correlation function in the grand-canonical ensem- 
ble then satisfies 

n s co 

(£ y Eo?o?)(n)Z(p) = J2 n s N XN Z N eM-»N). (4) 

a=l i,j N=l 

Using standard arguments [[HJ one can determine how the scaling behav- 
ior of the integrated spin-spin correlation function changes under coupling 
to gravity. In flat space, the spin-spin correlation functions scales as 

(of of) ~ \f i -f j \-^ A ' + ^\ (5) 
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The weight A° = A° = 1/16 is dressed by gravity in a theory of central 
charge c according to the KPZ formula 

(A CT - A°) = (l + 1 (VT^ - ) VT^j A ff (l - A CT ). (6) 

This dressed weight determines the scaling of the integrated spin-spin cor- 
relation function with fi on surfaces of genus h: 

EE^W%) - (a* - ao*- 1 ^™ 2 -^ 1 -*' (?) 



(the —1 before the dressed weight accounts for the integrations of i and j 
over the surface) with 



7s 12 



^(c-1- V(25-c)(l-c)). (8) 
The above relations and 

Z N ~ ;v- 1+ fr'- 2 X 1 - fc > (9) 
yield the finite-size scaling relation 

X* ~ JV^** (10) 

with 

7 /i/d H = 1 - 2A ff . (11) 

This is the scaling law that we shall verify numerically. The susceptibility 
scales as x ~ (/? — /5 C )~ 7 , the correlation length (governed by the decay of the 
spin-spin correlation function) obeys £ ~ (/3 — /3 C )~ 1 ' and dn is the intrinsic 
Hausdorff dimension of the random surface being considered. 

Assuming the standard scaling hyperscaling relation a = 2 — udn, we 
can also predict the value of 7. The specific heat scales as N a ' vdH . Then by 
applying the reasoning used to arrive at ([11]) to the two-point function of the 
energy operator e, one finds ajvdu equals 1 — 2A e , where A e is the dressed 
weight of the energy operator. It then follows that 

7= ( 1 - 2A -), (12) 
7 1-A e ) 1 > 
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One obtains A £ through the KPZ formula (B), substituting the bare energy 
weight A° = 1/2 for A°. 

We shall also measure scaling properties of the Fortuin-Kasteleyn (FK) 



clusters JT6| which we construct to update the spin degrees of freedom. The 



FK clusters appear in the reformulation of the Ising model as a correlated 



spin-bond percolation model with partition function [17 



z= E E p b a-pf b - b - (is) 

<T;=±1 colorings 

Colorings consist of a set of 'black' bonds drawn between adjacent points 
with identical spin values; each black bond is drawn with probability p = 
1 — exp(— 2/3). In (" p~3|) b bonds out of a possible N b bonds of the lattice are 
colored black. FK clusters comprise sets of sites linked together by black 
bonds. Therefore each cluster is assigned a single spin value. In the multi- 
generation case, we build a set of colored bonds and FK clusters separately for 
each species of spin. One can show that the spin-spin correlation function 
of the spins (crfcr") equals the pair-connectedness function (<5q\c Q ) of the 
corresponding FK clusters; 5 is 1 when i and j lie in the same cluster C a 



and otherwise |L8]|. From this, it follows that the spin-susceptibility equals 
the mean cluster size S = (s 2 )/(s), in which averages are taken over the 
distribution n(s), the mean number of clusters per configuration containing 
s sites. We thus shall determine the scaling of the mean cluster size and 
in addition, the fractal dimension df of the largest cluster. The average 
maximal size cluster M. of each configuration scales as 

M~N^. (14) 



Standard scaling arguments [HJ relate ^jvdu = 2df/dn — 1. One can derive 
this, for instance, by considering the asymptotic form of n(s) ~ Ns~ T . The 
standard hyperscaling relation vdu = 2 — a is then needed. The singularity of 
the cluster number density, the zeroth moment of n(s), as a function of p—p c , 
is given by the exponent (2— a). Similarly the second moment of n(s) scales as 
{p—Pc) 1 - Thus the usual scaling assumptions imply -y/udn = (r— 3)/(l— r). 
M. asymptotically obeys 

N 5- ~ 1 (15) 

J M 
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(that is, the mean number of clusters per configuration of size greater than 
M. is of order unity) and hence df /dn = l/(r — 1). Eliminating r then gives 
the above relation between df/dn and ^jvdu. 

We can also obtain additional information about the critical geometry 
of these theories by examining the properties of pure percolation clusters. 
Consider the bond-percolation model 

Z= E v\l-v) Nb ~ b q N % (16) 

colorings 

which for q = 2 is yet another formulation of the Ising partition function 
and more generally is the partition function of the g-state Potts model. The 
pair-connectedness function then exhibits the scaling behavior at criticality 



of (o) with weights [ 20 



*< = *. = W$y -(f) = (17) 



The q — > 1 limit corresponds to pure percolation, which has no dynamics 
(and vanishing central charge) and thus does not induce any back-reaction 
when it is coupled to a theory of gravity and matter of central charge c. 
A° j = 5/96; the dressed weight is again governed by the KPZ formula (f|). 
With this new dressed weight, we can then predict the scaling behavior of 
the mean (and maximal) cluster sizes SN,q=i (and M.N >q =i) using ([U]) with 
S substituted for x- 

In our simulations, we shall consider site percolation, in which sites 
(rather than bonds) are colored black with probability p and clusters are 
built by connecting adjacent colored sites. It is well known that bond and 
site percolation are in the same universality class, so that the scaling predic- 
tions described above should still hold in the case of pure site percolation. It 
is advantageous to consider site percolation because the critical value of p, 
p c , is constrained to equal 1/2 for site-percolation on triangulations PT|p[ 



2.1 c = 1 

The scaling relations become more complicated for c = 1. Analytic solutions 
of the c = 1 matrix models (and a careful analysis of Liouville theory) show 

2 There are some possible exceptions to this constraint, which definitely do not apply 
in the cases we shall consider here. This issue is discussed extensively in p4[. 
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that correlation functions no longer scale simply as powers of the cosmological 
constant //. Instead, the appropriate scaling variable is 77 which satisfies | 

fj, — — rj 01(77) + c\7] H (18) 

in the limit of small 77; c\ is a constant that we do not specify and shall 
not assume to be universal. This scaling relation has been derived for the 
Gaussian theory (of finite and infinite radius) coupled to gravity. The product 
of Ising models lies on the Gaussian c = 1 orbifold line |22] ; coupled to gravity, 



it is not equivalent to the solved matrix models. We shall assume that the 
asymptotic logarithmic scaling violation is characteristic of c = 1 and thus 
holds in the two-species case. Then we conjecture that the scaling relation 
(0) should be modified so that Q 

Q2o?of)Qi)ZQi) ~ ^M+a ct )+( 2 - 7s) (i-,)_ (19) 

In the case of pure percolation, the pair-connectedness function should be 
substituted for the spin-spin correlation function, and the dressed Ising 



weight A CT = 1/4 should be replaced by A CT)g=1 = ^5/96. In the follow- 
ing formulae, \ and the mean cluster size S are interchangeable. 

Our simulations will be done on worldsheets of toroidal topology (h = 1) 
for which Z{p) ~ ln(r/) for c = 1 (for those models that have been solved 
analytically, so again we are making an assumption about universality). To 
extract the asymptotic scaling behavior of xn, we invert the relation between 
77 and n order by order in 1/ ln/i and ln(— ln(/x))/ ln(/i) to obtain 

, = _^fi + MzM + f!2tM)'_H = M + ... 1 . (20) 

ln/i I ln/i \ m/x J (ln/i) 2 



We then expand the inverse Laplace transform of (|19f) to obtain 

I ( \n\nN f\n\nN\ 2 \n\nN 

V X vZv .V(.Vln.VF [ ' "hTA" • I "hTA" J " (LTiVT + " " ' 1 " 

/ u*(-u) r , AnlnN \ , M . 



3 Without loss of generality, we set fi c = 0. 

The essential role of logarithmic corrections in interpreting numerical measurements 
of 7 S at c = 1 has been previously discussed in reference p3j . 
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uj = 2(— 1 + A CT ) = —■y/udn — 1 (or = 2(— 1 + A (Tjg =i) for pure site 
percolation) and ^ is the digamma function. The scaling behavior of NZn 
is obtained by inverse Laplace transforming dZ{j]{^))/d^: 




lnlniV \ , . 

(i^lvF + '" ■ < 22 > 



In addition to the higher order terms that we have dropped from the inver- 
sion, there are additional corrections to the above formulae. The correction 
that depends on c\ in (|I5D, which we neglect, should lead to contributions 
to ([H]) and (^) that are competitive with the smallest corrections that we 
have included above. In addition, there should be the usual corrections to the 
scaling (|19D, but these will be suppressed by powers of 1/N and are negligible 
for moderately large N. 

The difference in the theoretically predicted scaling behavior of percola- 
tion clusters at c = 1, compared to c = 1/2, illustrates the sensitivity of the 
worldsheet geometry to the presence of the Ising spins. This should induce an 
effective coupling between species, which we should be able to detect through 
correlations between different species' observables. To look for this coupling, 
we shall measure 

e a p = ±((ErEfl)-(E^(Efi)) (23) 

and 

m a p = ±({\M<*\\MP\)-{\M<*\){\MP\)) (24) 

with E a = Cij (r)crfcrj and \M a \ — J2i \°f\- A nonzero effective coupling 
between species should then be manifest through the quantities e* = 2e 12 /tr 
e and m* = 2m 12 /tr m, which vanish when E a (\M a \) are uncorrelated and 
are 1 or —1 when they are respectively perfectly correlated or anti-correlated. 



Numerical Simulations and Results 



The partition function ([[]) is evaluated numerically by a Monte Carlo simu- 
lation. The sum over triangulations is implemented via the standard DTRS 
algorithm [53|. This updates the connectivity matrix CV,- by flips on ran- 
domly chosen pairs of triangles sharing a common link. The Ising spins are 
updated using the Swendsen-Wang algorithm P5|; one builds FK clusters 
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over the entire lattice and then assigns a randomly chosen value to the spin 
associated with each cluster. 

Runs were performed on lattices of toroidal topology of size N = 2048, 4096, 8192 
and 16384. Each sweep consisted of 3iV flips of randomly chosen links fol- 
lowed by a Swendsen-Wang update of the spins. To perform the percolation 
measurements the sites on the lattice were then randomly colored and perco- 
lation clusters constructed. We used the jacknife technique to estimate our 
errors. We measured auto-correlation functions and computed integrated 
auto- correlation times, using standard techniques |i~7| . The magnetization, 
susceptibility and cluster sizes exhibited considerable critical slowing down 
0. It was necessary therefore to sample quite a large number of lattices 
- for each data point between 15, 000 and 30, 000 independent lattices, re- 
quiring from 30, 000 to 900, 000 sweeps. The total CPU-time used was ap- 
proximately equivalent to six months on an HP-9000 (720) workstation. We 



histogrammed our data [^7]]. For cluster data on the larger lattices, how- 
ever, histogramming was not reliable, given our statistics. The cluster data 
exhibited extremely large fluctuations from one measurement to the next; 
presumably this was the source of the poor performance of histogramming. 

To extract critical exponents using finite-size scaling we must first deter- 
mine the critical point. We estimated the value of f3 c by locating the peak in 
the lattice susceptibility (3/N((M 2 ) - (|M|) 2 ), where M is the magnetization 
averaged over species []. We first present our results for a single species. Fig- 
ure [l] reveals that this quantity peaks at a value of (3 C = .2185(20). We also 
plotted Binder's cumulants @ U M = 1 - (M 4 )/(3(M 2 ) 2 ) as a function of (3 
for different lattice sizes. From the position of the intersections of these cu- 
mulants, we estimate (3 C = .2180(7). The critical temperature for Ising spins 
on lattices dual to our triangulations was computed analytically in [p9fl , so 
by the Ising duality relation, it follows that /3 C = (1/2) ln(131/85) ~ .216273. 
Our numerical estimate of the critical value of (3 is therefore somewhat high. 

One could also envision locating the critical temperature by looking for a 
minimum of the mean pure percolation cluster size as a function of (3. The 
growth in the mean cluster size is determined by the exponent (7/1/^)^=1, 

5 More detail on critical slowing down for various algorithms applied to spin models on 
random lattices will appear in future work p6| . 

6 Note that this does not equal the susceptibility xn defined in section (||). xn contains 
no subtractions and agrees with the continuum susceptibility only in the high-temperature 
phase. 
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Figure 1: The 'lattice susceptibility' plotted for the one-species model for 
lattices of size 1024 — 16384. Error bars are only shown when they are larger 
than the symbols. 
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0.215 



0.216 



0.217 



0.218 



0.219 



Figure 2: Binder's cumulants Um for N = 2048,4096 and 8192 lattices in 
simulations of the one species Ising model for lattices of size 1024 — 16384. 
Lines are drawn to guide the eye. 
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320 1 1 1 ' 1 1 1 1 1 1 1 1 1 ' 

0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 

P 



Figure 3: The mean size of percolation clusters as a function of (5 for the 
one species Ising model on 2048 node lattices. The dip occurs near criticality. 
The dotted line is drawn to guide the eye. 

which decreases with c. We would therefore expect that percolation clusters 
for a given lattice size will become smaller as c increases or likewise as (3 is 
tuned to bring the Ising spins to criticality. In Fig. we do indeed observe 
a dip in the mean cluster size around the estimated value of [3 C . The dip 
is quite broad, however, and doesn't pinpoint (3 C ; the distinction between 
the j3 = .218 and .210 points in the figure may be a statistical artifact. 
For N = 4096, for instance, with better statistics we find that the mean 
size is measured to be 571.8 ± 1.9,571.1 ± 1.9,573.1 ± 1.9,572.0 ± 1.8 and 
572.6 ± 1.9 for (3 = .211, .213, .215, .216273 and the .218. Thus we observe no 
significant variation over a range 5/3 = .007. On the other hand the breadth 
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N 


1024 


2048 


4096 


8192 


theory 


(~l/vd H ) eff from \ 


.702 (6) 


.705 (5) 


.686 (6) 


.688 (10) 


2/3 


(j/ud H ) ef f from S FK 




.701 (7) 


.690 (9) 


.676 (12) 


2/3 


{d f /d H ) eff 




.844 (5) 


.839 (7) 


.830 (8) 


5/6 






.747 (6) 


.741 (7) 


.722 (10) 


.710 


{df/d H )% 




.866 (5) 


.866 (6) 


.852 (8) 


.855 



Table 1: Summary of exponents extracted from finite-size scaling at the exact 
value of P c for one Ising species. The exponents are computed using data 
from lattices of size N and 2N. 



of the dip means that measurements of / vdn) q= i are not very sensitive 
to the estimated value of j3 c . This is not true for the exponent (7/1/^) 
which, as we shall see, exhibits a strong dependence on the estimated critical 
temperature. 

From the finite-size scaling relation (|10|) we define the effective exponents 

( T /^) e// =lnf^/ln2. (25) 
V Xn ) 

When cluster sizes rather than the spin-susceptibility are considered, it is 
implicit that S is to be substituted for x m the above definition. Likewise, 



(d f /d H ) eff = 




As usual, the analogous exponents for pure percolation are defined as above 
with subscript q — 1. We now summarize our results for these critical expo- 
nents, extracted through finite size scaling of xn, the mean FK cluster size 
Sfk, the maximal FK cluster size Mfk, the mean percolation cluster size 
S q =i and the maximal percolation cluster size M. q= \. In table ([[]), we present 
these values obtained from runs at the known value of (3 C ~ .216273. The 
exact values are "i/vd H = 2/3, df/d H = 5/6, (7/z/d^) g=1 = (4— y7/2)/3 and 

(df /dH) q =i = (7— ^7/2)/6. The agreement between our measurements and 
the theoretical predictions is quite good. Since (df/dn) only describes prop- 
erties of the largest cluster, it might be less subject to corrections to scaling 
(for small cluster sizes) and thus provide the best estimator of the critical 
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N 


1024 


2048 


4096 


8192 


theory 


(~l/vd H ) eff from \ 


.74 (3) 


.76 (3) 


.73 (2) 


.76 (4) 


2/3 


(j/ud H ) ef f from S FK 




.748 (14) 


.741 (8) 




2/3 


(d f /d H ) eff 




.877 (15) 


.869 (6) 




5/6 


{l/vd)eff,q=l 




.750 (10) 


.739 (8) 




.710 


(df/d H )eff,q=l 




.871 (7) 


.867 (5) 




.855 



Table 2: Summary of exponents extracted from finite-size scaling using the 
numerically determined value of (3 C for one Ising species. 



exponents. Indeed, the values of (df/d^eff for FK and percolation clus- 
ters on the larger lattices already match the asymptotic predictions, within 
our statistics. The exponents (7 jvdu)efj decrease towards their asymptotic 
values and differ from them by only about 2% on the largest lattices. 

In the two species case we do not, however, know (3 C . To use the single 
species model as a control, therefore, we estimate the critical exponents using 
finite-size scaling at an estimated value of f3 c = .2180(7). We present these 
estimates in table [2|. The values quoted from the scaling of x characterize 
the range of exponents in the window (3 G [.2173, .2187]; this determination 
relies on histogramming. The exponents extracted from cluster data are 
taken from data at (3 = .218 (within this window), since histogramming was 
not reliable for this data. We see that the shift in temperature away from 
the exact (3 C induces a large change in the Ising exponents. In particular, 
[^ljvdu)eff^ which previously came within 2% of the asymptotic value on 
the largest lattice, now is 10-15% higher. The percolation exponents fare 
much better; they are quite insensitive to the value of (3 and agree with the 
values taken at j3 c = .216273. 

We also attempted to perform a direct fit of x to (/? — /5 C )~ 7 . We found 
that a power law fit only seemed to work in the region /3 = .19 — .195, where 
finite size effects were not severe. Only two of our data points (corresponding 
to the lowest values of (3 at which we ran) were then used in this fit, yielding 
7 = 1.8(1), which is not so far from the theoretical value of 7 = 2 that 
follows from ([12]). Given the sparsity of our data in this regime and thus 
the difficulty in verifying that we are seeing asymptotic scaling, these results 
should be interpreted with caution. 

We now turn to a treatment of the two species results. First, we present 
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Figure 4: A plot of the 'lattice susceptibility' for the two species case with 
N = 2048,4096,8192 and 16384. Error bars are only shown when they are 
larger than the symbols. 

plots of P/N({M 2 ) — (|M|) 2 ) and the intersections of Binder's cumulants. 

Note that these observables look qualitatively very similar to their one 
species counterparts. The position of the susceptibility peak is (5 = .2185^20- 
The Binder's cumulants (taken from the largest lattices where we have suffi- 
cient data to determine intersections) give f3 c = .217(1). As in the one-species 
case we adopt the Binder's estimate, which seems to be more precise. Based 
on the logarithmic corrections to scaling discussed in section [| we would es- 
timate that on the lattice sizes simulated we are further from the asymptotic 
scaling regime than in the one species case. It is therefore likely that our 
estimate of (3 C here is less precise. 

In table H, we compare the finite-size scaling measurements of (j/vd^eff 
and (7 / udnYgdi with the theoretical exponents, computed from the KPZ for- 
mula, neglecting logarithmic corrections. The first row contains an estimate 
of the susceptibility exponent based on histogram extrapolations within the 
region (3 = .216 to .218. We observed that measuring the mean-size of FK 
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0.61 



0.605 - 



0.595 



0.59 



0.585 




0.58 - 



0.575 



0.215 



0.216 



0.217 



0.218 



0.219 



Figure 5: Binder's cumulants for N = 4096 and N = 8192 as measured in 
the two species model. 





1024 


2048 


4096 


8192 


theory 


(l/"d)eff from x 


.73 (3) 


.73 (3) 


.75 (5) 


.74 (4) 


1/2 


(l/vd) efLq=1 


.759 (11) 


.717 (5) 


.712 (7) 


.709 (10) 


.544 


(df/d H )eff,q=l 


.852 (6) 


.848 (4) 


.848 (5) 


.848 (6) 


.772 


(7/vd) efLq=1 (X) 


.739 (7) 


.732 (7) 


.700 (10) 




.544 


(df/ d H ) e ff s= i(X) 


.861 (5) 


.859 (7) 


.841 (8) 




.772 



Table 3: Summary of exponents extracted from finite-size scaling for two Ising 
species compared with KPZ exponents, without logarithmic corrections. In 
the last two rows, we present the corresponding percolation exponents for 
lattices on which a Gaussian field X was simulated. 
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clusters (at (3 = .216) always yielded exponents agreeing with those extracted 
from the spin-susceptibility. Since the cluster data did not histogram reliably, 
it was difficult though to estimate the corresponding exponents throughout 
the above range of (3 without taking a large amount of additional data. FK 
cluster scaling was primarily measured just to verify that it agreed with the 
scaling of the spin-susceptibility. Our data already showed this, so it did 
not seem worthwhile to repeat runs to measure the cluster scaling at various 
temperatures. The following rows thus summarize data for pure percolation 
clusters, which was relatively insensitive to (3; i.e. a shift of (3 of .005 did 
not induce a stastically significant change in the q = 1 exponents. Since 
again percolation cluster data did not histogram reliably, we present data 
taken at (3 — .216 in table |3|; these should be representative of the values one 
would measure throughout the region (3 G [.216, .218]. The final two rows 
include exponents extracted from simulations of a Gaussian field X coupled 
to gravity, as discussed in [ITJ . 



What is most striking about the measured c = 1 exponents is in fact 
that they agree quite well with the c = 1/2 data (taken from the numerically 
estimated range of (3 C ) and in the case of percolation, the c = 1/2 theoretical 
predictions. There is clearly, however, a very large discrepancy between the 
measured and theoretical c = 1 scaling exponents. On the larger lattices, the 
percolation exponents for the two species Ising model also agree fairly well 
with those measured in the Gaussian field simulations, suggesting some de- 
gree of universality at c = 1 in their behavior. From equations fl2jj) and (|22]) , 
we can determine the logarithmic corrections to the exponents (7 / 'vdjjjeff 
and inM H )% One can easily see that these are considerable by including 
the leading correction, which gives 

(TMr).// = i/vd B + 1 + ^ dH +■■■; (27) 

this formula of course holds for q — 1. In Fig. ^| we compare the measured 
{'y/vdnjeff with the corresponding theoretical predictions including both the 
leading correction ( P7|) and all of the logarithmic corrections that follow from 
equations fl2T|) and (|22|). The data does not agree particularly well with any of 
the predictions, but that is to be expected, since our estimate of f3 c most likely 
induces a large error. Note that the theoretical curve including logs lies very 
close to the c = 1/2 asymptotic value of j/vdu = 2/3 for the lattice sizes we 
simulate. The measured values of this exponent for c = 1 and c = 1/2 (when 
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0.85 




2000 4000 6000 8000 

Lattice size N for N to 2N extrapolation 



Figure 6: A comparison of the theoretical prediction of the finite-size c = 
1 magnetic susceptibility scaling with our data. The dotted line includes 
the leading logarithmic correction and the solid line takes into account all 
subleading terms that we calculated. The horizontal dashed line represents 
the prediction without logarithmic corrections. 
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we use the numerically estimated values of (3 C ) are essentially identical, so 
that the discrepancy between the theoretical predictions (incorporating logs 
in the c = 1 case) and the data are thus roughly the same for c = 1 and 
c= 1/2. 

In Fig. [7] a similar comparison is shown for the exponent (7 jvdn) e ^l x . We 
see that the data and theoretical predictions match quite well. Presumably, 
as in the c = 1/2 case, the comparison works because this exponent no longer 
depends sensitively on our determination of (3 C . The agreement with theory is 
about as successful as in the Gaussian case Rl4| . This suggests the likelihood 
that at least the leading logarithmic corrections to scaling (^) are correct 
and universal for c = 1. We should note here that although the data can be 
fit by including the above logarithmic corrections it would not be possible to 
extract these corrections from the data itself without theoretical guidance. 

As in the one species case, we also fitted our lowest (3 data points to 
X ~ (P — /3 C )~ 7 - The fit yielded 7 = 2.03(4) which does not match the 
theoretical value of 7 = 1 + 1/V2 ~ 1.71. It is not clear that this fit is 
reliable; probably the exact relation between x an d P should also include 
logarithmic corrections. 

We finally turn to an examination of correlations (defined in and 
fl23"D ) between the different spin species. Fig. || exhibits definite, though 
moderately small, correlations between the spins of different species. As 
evident in Fig. the correlations between the average species' energies, as 
measured by e*, are much smaller. The magnetization (and more directly, 
the susceptibility) is strongly correlated with the distribution of FK cluster 
sizes, which in turn should be sensitive to the bottlenecks which characterize 
the worldsheet geometry. Since the correlations between species are mediated 
by fluctuations in the geometry, it is not surprising that they are stronger in 
the magnetic sector than in the energy sector. In both cases the correlations 
are not so strong, indicating that the Ising spins are only weakly coupling to 
gravity for the lattice sizes we consider. 



4 Discussion 

To gauge our prospects for extending this work to c > 1, we shall now 
attempt to shed more light on the origin of the corrections to scaling at 
c = 1. The dynamics at c = 1 is governed by the tachyon (the lowest mass 
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Figure 7: A comparison of the theoretical prediction of the finite size c = 1 
scaling of (7/W.ff)g=i with our percolation data for the two species Ising 
model. The dotted line includes the leading logarithmic correction and the 
solid line takes into account all subleading terms that we calculated. The 
horizontal dashed line indicates where these curves asymptote. 
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0.222 



Figure 8: The correlation m* plotted for the two species model on N = 8192 
lattices. Recall that the correlations are normalized so that m* equals one 
when |M Q=1 | and |M a=2 | are fully correlated. 
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Figure 9: The correlation e* between the energies of different species plotted 
for the c = 1 model on N = 8192 lattices. 



operator) which acquires a continuous spectrum with vanishing energy at 
momentum. The scaling corrections arise from the infrared cutoff on the 



tachyonic momentum, p > 1/ In AT [ . The origin of this logarithmic cutoff 
may not be so evident in the context of the two-species Ising model, since 
the tachyon cannot be expressed so simply in terms of the Ising spins. The 
presence of logarithmic corrections appears to only depend on the value of 
the central charge, so we appeal to universality and consider instead the c = 1 
Gaussian theory. The tachyon can then be written in terms of the Gaussian 

field X as 

T{p) ~ / M ex V (ipX + (\p\ - 2)0); (28) 



is the Liouville field and g is the reference metric (g = gexp(— 0)). The 
infrared target space momentum cutoff is then just a consequence of the 
infinite Hausdorff dimension of the embedding space p0| , 



(XX) ~ (In Nf 



(29) 



The suppression of tachyon propagation by finite-size effects effectively weak- 
ens the coupling between gravity and the Ising spins. Therefore, the spins 
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of each species are only very weakly coupled, and the two and one species 
models appear to be qualitatively very similar on the lattices we consider. 
This should not be true in the continuum limit. We see that without an 
understanding of the corrections to scaling for c = 1, the numerical observa- 
tions seem to conflict with our theoretical expectations of the onset of strong 
coupling at c = 1. 

The target space Hausdorff dimension, as measured for Gaussian embed- 
dings of d somewhat greater than 1, is also large. Thus we would expect 
that for c > 1, the tachyon ground state energy still acquires a considerable 
shift and we expect large corrections to scaling. At c = 1, we were able 
to determine these corrections because the c = 1 model is solvable. With- 
out new theoretical input, the prognosis for understanding simulations for c 
somewhat greater than 1 thus seems poor. 
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